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Abstract 



Prolate spheroidal wave functions (PSWFs) play an important role in various areas, from 
physics (e.g. wave phenomena, fluid dynamics) to engineering (e.g. signal processing, filter de- 
sign). One of the principal reasons for the importance of PSWFs is that they are a natural and 
efficient tool for computing with bandlimited functions, that frequently occur in the abovemen- 
tioned areas. This is due to the fact that PSWFs are the eigenfunctions of the integral operator, 
that represents timelimiting followed by lowpassing. Needless to say, the behavior of this op- 
erator is governed by the decay rate of its eigenvalues. Therefore, investigation of this decay 
rate plays a crucial role in the related theory and applications - for example, in construction of 
quadratures, interpolation, filter design, etc. 

The significance of PSWFs and, in particular, of the decay rate of the eigenvalues of the 
associated integral operator, was realized at least half a century ago. Nevertheless, perhaps 
surprisingly, despite vast numerical experience and existence of several asymptotic expansions, 
a non-trivial explicit upper bound on the magnitude of the eigenvalues has been missing for 
decades. 

The principal goal of this paper is to close this gap in the theory of PSWFs. We analyze the 
integral operator associated with PSWFs, to derive fairly tight non-asymptotic upper bounds on 
the magnitude of its eigenvalues. Our results are illustrated via several numerical experiments. 

Keywords: bandlimited functions, prolate spheroidal wave functions, eigenvalues 
Math subject classification: 33E10, 34L15, 35S30, 42C10, 45C05, 54P05 

1 Introduction 

The principal purpose of this paper is to establish and prove several inequalities involving the 
eigenvalues of a certain integral operator associated with bandlimited functions (see Section|3]below). 
While some of these inequalities are known from "numerical experience" (see, for example, [4], [9], 
|15|). their proofs appear to be absent in the literature. 

A function / : R — > M. is bandlimited of band limit c > 0, if there exists a function a € L 2 [—1, 1] 
such that 



In other words, the Fourier transform of a bandlimited function is compactly supported. While ([T]) 
defines / for all real x, one is often interested in bandlimited functions, whose argument is confined 

*This author's research was supported in part by the AFOSR grant #FA9550-09-l-0241 
TYale University, 51 Prospect st, New Haven, CT 06511. Email: andrei.osipov@yale.edu. 




(1) 



1 



to an interval, e.g. — 1 < x < 1. Such functions are encountered in physics (wave phenomena, fluid 
dynamics), engineering (signal processing), etc. (see e.g. [T3], [IS], [2"0]). 

About 50 years ago it was observed that the eigenfunctions of the integral operator F c : L 2 [— 1, 1] — > 
L 2 [—1,1], defined via the formula 



provide a natural tool for dealing with bandlimited functions, defined on the interval [—1,1]. More- 
over, it was observed (see 0, 0, [TT]) that the eigenfunctions of F c are precisely the prolate 
spheroidal wave functions (PSWFs), well known from the mathematical physics (see, for exam- 
ple, [TB], [!!])■ The PSWFs are the eigenfunctions of the differential operator L c , defined via the 
formula 



In other words, the integral operator F c commutes with the differential operator L c (see 0, |18|). 
This property, being remarkable by itself, also plays an important role in both the analysis of PSWFs 
and the associated numerical algorithms (see, for example, 0, [3]). 

Obviously, the behavior of the operator F c is governed by the decay rate of its eigenvalues. Over 
the last half a century, several related asymptotic expansions, as well as results of numerous numerical 
experiments, have been published; moreover, implications of the decay rate of the eigenvalues to both 
theory and applications have been extensively covered in the literature - see, for example, 0, 0, 
H- 0, 0, 0,0, US], QI], P2> P3]> B3> B3- Tt is perhaps surprising, however, that a non-trivial 
explicit upper bound on the magnitude of the eigenvalues of F c has been missing for decades. This 
paper closes this gap in the theory of PSWFs. 

This paper is mostly devoted to the analysis of the integral operator F c , defined via (J2J). More 
specifically, several explicit upper bounds for the magnitude of the eigenvalues of F c are derived. 
These bounds turn out to be fairly tight. The analysis is illustrated through several numerical 
experiments. 

Some of the results of this paper are based on the recent analysis of the differential operator 
L c , defined via (0), that appears in [22], [23]. Nevertheless, the techniques used in this paper are 
quite different from those of [22] , [23] . The implications of the recent analysis of both L c and F c to 
numerical algorithms involving PSWFs are being currently investigated. 

This paper is organized as follows. In Section [2] we summarize a number of well known math- 
ematical facts to be used in the rest of this paper. In Section [3j we provide a summary of the 
principal results of this paper, and discuss several consequences of these results. In Section 2] we 
introduce the necessary analytical apparatus and carry out the analysis. In Section [5] we illustrate 
the analysis via several numerical examples. 

2 Mathematical and Numerical Preliminaries 

In this section, we introduce notation and summarize several facts to be used in the rest of the 
paper. 

2.1 Prolate Spheroidal Wave Functions 

In this subsection, we summarize several facts about the PSWFs. Unless stated otherwise, all of 
these facts can be found in 0, 0, 0, 0, 0, [22], [23J. 
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Given a real number c > 0, we define the operator F c : L 2 [— 1, 1] — > L 2 [— 1, 1] via the formula 

F c [<p] (x) = J ^(t)e lcxt dt. (4) 

Obviously, F c is compact. We denote its eigenvalues by Ao, Ai, . . . , A„, . . . and assume that they 
are ordered such that |A n | > |A„+i| for all natural n > 0. We denote by tp n the eigenfunction 
corresponding to A n . In other words, the following identity holds for all integer n > and all real 
-1 < x < 1: 

Ki>n(x) = J Mty^dt. (5) 

We adopt the convention^ that |] V^n 1 1 i,2[— i, i] = 1- The following theorem describes the eigenvalues 
and eigenfunctions of F c (see [3J, [I], [H])- 

Theorem 1. Suppose that c > is a real number, and that the operator F c is defined via (|4]) 
above. Then, the eigenfunctions tpo^tpx, ... of F c are purely real, are orthonormal and are complete 
in L 2 [—1, 1]. The even-numbered functions are even, the odd-numbered ones are odd. Each function 
ip n has exactly n simple roots in (—1,1). All eigenvalues X n of F c are non-zero and simple; the 
even-numbered ones are purely real and the odd-numbered ones are purely imaginary; in particular, 
A„ = i"|A n |. 

We define the self-adjoint operator Q c : L 2 [— 1, 1] — > L 2 [— 1, 1] via the formula 

p -. . . If 1 sm(c(x — t)) , . , , . 

Q c [<p] (x) = - ' <P(t) dt. (6) 

7T J_ x X-t 

Clearly, if we denote by 2f : L 2 (R) — > L 2 (M.) the unitary Fourier transform, then 

QcM(x) = X[ -i,i](x)-?- 1 [x[-c,c){0-?[v}(0] (x), (7) 
where X[~a,a] '■ K — >■ M is the characteristic function of the interval [—a, a], defined via the formula 

( ) < ^ ~ a — x — a ' ($\ 
' a ' a ' 1 otherwise, 

for all real x. In other words, Q c represents low-passing followed by time-limiting. Q c relates to F c , 
defined via (gj), by 

Qc = ■ F* ■ F c , (9) 

Z7T 

and the eigenvalues fi n of Q n satisfy the identity 

M„ = ^-|A„| 2 , (10) 
for all integer n > 0. Moreover, Q c has the same eigenfunctions ip n as F c . In other words, 

, , . 1 f 1 sin(c(x — t)) , , . , 

ihcK 0* = - / — , Mt) dt, (n) 

7T X-t 



1 This convention agrees with that of [3] , 0] and differs from that of [5] . 



3 



for all integer n > and all — 1 < x < 1. Also, Q c is closely related to the operator P c : i 2 (R) — > 
L 2 (R), defined via the formula 

„ r i / x 1 /"°° sin (c(x — t)) , . . , „. 

P c ^ x = - / V V jj ¥»(t) dt, 12 

which is a widely known orthogonal projection onto the space of functions of band limit c > on 
the real line K. 

The following theorem about the eigenvalues [i n of the operator Q Cl defined via ©, can be traced 
back to 

Theorem 2. Suppose that c > and < a < 1 are positive real numbers, and that the operator 
Q c : L 2 [— 1,1] — > [— 1,1] is defined via © above. Suppose also that the integer N(c,a) is the 
number of the eigenvalues fi n of Q c that are greater than a. In other words, 

N(c, a) = max{fe = 1, 2, . . . : fik-i > a } ■ (13) 

Then, 

N(c,a) = — + (%logi^V°g c + °( lo g c )- ( 14 ) 
7r \ir z a J 



According to (|14l) . there are about 2c/7r eigenvalues whose absolute value is close to one, order of 
log c eigenvalues that decay exponentially, and the rest of them are very close to zero. 

The eigenfunctions ip n of Q c turn out to be the PSWFs, well known from classical mathematical 
physics [TB]. The following theorem, proved in a more general form in [TT], formalizes this statement. 

Theorem 3. For any c > 0, there exists a strictly increasing unbounded sequence of positive numbers 
Xo < Xi < • • • such that, for each integer n > 0, the differential equation 

(1 - x 2 ) ■ ip"(x) - 2x ■ ip'(x) + (x-n - c 2 x 2 ) ■ ip(x) = (15) 

has a solution that is continuous on [—1,1]. Moreover, all such solutions are constant multiples of 
the eigenfunction ip n of F c , defined via ^) above. 

In the following theorem, that appears in [4], an upper bound on |A n | in terms of n and c is 
described (the accuracy of this bound is discussed in Section 13.21 below; see also Theorem [34] and 
Remark [TT1 in Section [ 



Theorem 4. Suppose that c > is a real number, and n > is a non-negative integer. Suppose 

also that X n is the nth eigenvalue of the operator F c , defined via Suppose furthermore that the 
real number v(n, c) is defined via the formula 

0F-c" (n!) 2 

Kn ' c)= (2n)M> + 3/2)' (16) 
where T denotes the gamma function. Then, 

|A„| <i/(n,c). (17) 

Moreover, 

X n (c) = i n u(n,c)-e R ^ c \ (18) 
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where the real number R(n, c) is defined via the formula 



R(n, c) 



2(^(1)^-1 n 
2t t 



dr. 



(19) 



The function ip^ in (|19[) is the nth PSWF corresponding to the band limit r. 

The following approximation formula for |A n | appears in Theorem 18 of d], without proof (though 
the authors do illustrate its accuracy via several numerical examples). 

Theorem 5. Suppose that c > 1 is a real number, and that n > c is a positive integer. Suppose also 
that the real number po(n, c) is defined via the formula 



p (n,c) 



■ exp 



fa' F 



•Xr, 



' Xn ~ C z 



Xr, 



Xn 



(20) 



where F, E are the complete elliptic integrals, defined, respectively, via (|39| , (|40| in Section \2.3[ 
Then. 



\Xt. 



Po{n,c) 



O 



1 



(21) 



Remark 1. Obviously, (121|) cannot be used in rigorous analysis, due to the lack of both error esti- 
mates and proof. In addition, the assumption n > c turns out to be rather restrictive. Nevertheless, 
in Section^we establish several upper bounds on \X n \, whose form is similar to that of po(n,c). 
The approximate formula (|2ip will only be used in the discussion of the accuracy of these bounds, 
in SectionWM. 

The following four theorems contain relatively recent results. All of them appear in [22 , [23 . 

Many properties of the PSWF ip n depend on whether the eigenvalue Xn of the ODE (fT5|) is greater 
than or less than c 2 . In the following theorem from |22j . [23], we describe a simple relationship 
between c, n and \ n . 

Theorem 6. Suppose that n > 2 is a non-negative integer. 
•Ifn< (2c/ir) - 1, then \ n < c 2 . 

• If n ^ (2c/7r), then \n > c 2 - 

• If (2c/ n) — 1 < n < (2c/ir), then either inequality is possible. 



In the following theorem from |22) . |23) . we describe upper and lower bounds on Xn m terms of n 
and c. 

Theorem 7. Suppose that n > 2 is a positive integer, and that \n > c 2 . Then, 



7T 



where the function E : [0, 1] — > R is defined via (|40[) in Section 



1 fa'E[-^= ) < n + 3, (22) 
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In the following theorem, we provide another upper bound on Xn in terms of n. 
Theorem 8. Suppose that n> 2 is a positive integer, and that Xn > c 2 . Then, 

Xn< (| (n + 1)) 2 . (23) 

In the following theorem, we describe an upper bound on the reciprocal of |^ n (0)| for even n (see 
Theorem 21 in [22]). 

Theorem 9. Suppose that n > is an even integer, and that Xn > c 2 . Then, 



' <4-Jn-§. (24) 



l^(0)| " V c 

Remark 2. Detailed numerical experiments, conducted by the author, seem to indicate that, in fact, 

' 0(1) (25) 



l^n(0)| 

(see also RV). In other words, the inequality (|24j) is rather crude; on the other hands, it has been 
rigorously proved, and is sufficient for our purposes. 

2.2 Legendre Polynomials and PSWFs 

In this subsection, we list several well known facts about Legendre polynomials and the relationship 
between Legendre polynomials and PSWFs. All of these facts can be found, for example, in [7J, [3] 

EH- 

The Legendre polynomials Pq, Py, P%, . . . are defined via the formulae 

Po(t) - 1, 

P l {t)=t, (26) 

and the recurrence relation 

(k+l)P k+1 (t) = (2k + l)tP k (t)-kP k ^(t), (27) 

for all k = 1, 2, . . . . The Legendre polynomials {-Pfej^Lo constitute a complete orthogonal system in 
L 2 [—1, 1]. The normalized Legendre polynomials are defined via the formula 



Ph(t) = Pt(t) ■ Vfc+V2. (28) 

for all k = 0, 1, 2, The L 2 [—1, l]-norm of each normalized Legendre polynomial equals to one, 

i.e. 

1 

(K(t)) 2 dt = l. (29) 

Therefore, the normalized Legendre polynomials constitute an orthonormal basis for L 2 [— 1 , 1] . In 
particular, for every real c > and every integer n > 0, the prolate spheroidal wave function ip n , 
corresponding to the band limit c, can be expanded into the series 

oo 

i, n (x)=Y J ^' C) ■P'k(x), (30) 

fe=0 



G 



for all — 1 < x < 1, where /3o™' c \ Pi' c \ • • ■ are defined via the formula 



/3["' c) = J 4> n {x) ■ Pk(x) dx, (31) 



for all fc = 0, 1, 2, ... . The sequence /?o™' C \ /3[ n ' c \ . . . satisfies the recurrence relation 







= Xn 






^(n,c) 


Xn 




Afc,fc- 2 • 2 C) + A Kk ■ (3 { k n - c) + A kik+2 


o(n,c) 
' ^k+2 


= Xn 





(32) 

for all fc = 2, 3, . . . , where A k +2,k, Ak^k+2 arc defined via the formulae 

2fc(fc+l)-l 2 
^- fc ( fc + 1 ) + (2 fc + 3) (2 i-l) - C ' 
A -A (k + 2)(k + l) 

Ak.k+2 — Au+2,k — : — ■ C , (63) 

(2k + 3) v / {2k + l){2k + 5) 

for all fc = 0, 1, 2, In other words, the infinite vector f3 — |/3[™' c ' > | satisfies the identity 

(A - Xn I) ■ P = 0, (34) 
where the non-zero entries of the infinite symmetric matrix A are given via (I33[) . 

2.3 Elliptic Integrals 

In this subsection, we summarize several facts about elliptic integrals. These facts can be found, for 
example, in section 8.1 in [7], and in |21j . 

The incomplete elliptic integrals of the first and second kind are defined, respectively, by the 
formulae 

F(y,k)= / - (35) 
Jo VI - k 2 sin 2 t 

E(y, fc) = [ V Vl- fc 2 sin 2 1 dt, (36) 
Jo 

where < y < n/2 and < k < 1. By performing the substitution x — sini, we can write (|35[) and 

as 

F(y ' fc) y v^ 1 - z 2 ) I 1 -"^T (37) 

rsinfo) ll_ k 2 x 2 
l-X 2 

The complete elliptic integrals of the first and second kind arc defined, respectively, by the formulae 

F(k) = F (* fc) = - f (39) 

v 2 y Vo V 1 - k 2 sin 2 * 

/7T \ r /2 , 

E(k) = E^-,k) = j V 1 - fc 2 sin 2 t dt, (40) 



E (y> k )= I \l- 2-dx. (38) 
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for all < k < 1. Moreover, 

E (VT~k^) = 1 + (-1 + log(2) - • fc 2 + O (fc 4 • log(fc)) . (41) 

In addition, 

F{k) - E(k) > j ■ fc 2 , (42) 

for all real < k < 1. 

3 Summary and Discussion 

In this section, we summarize some of the properties of prolate spheroidal wave functions and the 
associated eigenvalues, proved in Section |U In particular, we present several upper bounds on 
|A„| and discuss their accuracy. The PSWFs and related notions were introduced in Section |2~T1 
Throughout this section, the band limit c > is assumed to be a positive real number. 



3.1 Summary of Analysis 

In the following two propositions, we provide some upper bounds on the eigenvalues \n of the ODE 
(fTS")) . They are proved in Theorem I2H HH1 m Section [ 



Proposition 1. Suppose that n is a positive integer, and that 



7T 7T 

for some 



2c 2 . , {4ewc\ 
n > — + -2 ■ S ■ log — — ) , (43) 



< 6 < ^ • c. (44) 



Then, 

4 



Xn > ? + - • 6 ■ c. (45) 

7T 



Proposition 2. Suppose that n is a positive integer, and that 



for some 



The 



2c 2c 2 p , /4evrc\ 

< n < - + ■ 8 ■ log — - 3, (46) 



3<5<--c. (47) 



2 8-6 

Xn<c 2 + c. (48) 

7T 
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The following is one of the principal results of this paper. It is proved in Theorem[33]in Section l4~2l 
(see also Remark [5]), and is illustrated in Experiments 2, 3 in Section [5j 

Proposition 3. Suppose that n > is an even integer number, and that A n is the nth eigenvalue 
of the integral operator F C) defined via ([5]) in Section \2.1\ Suppose also that 



2C /77T 

n > h V42. 

Suppose furthermore that the real number £(n, c) is defined via the formula 



C(n,c) 



7 



(4 • Xn/C 



2|^„(0)| 3- X n/c 2 -l 

(f 



(Xn - C 



2\ 4 



exp 



X, 



E 



Xr, 



Xr, 



(49) 



(50) 



where Xn is the nth eigenvalue of the differential operator L c , defined via ([3]) in Section]^ and F,E 
are the complete elliptic integrals, defined, respectively, via (|39p . (|40[) in Section ] 2. 3[ Then, 



| A„ | < ((n,c). 



(51) 



Remark 3. It follows from the combination of Remark^ in Section \2.1\ and Proposition^ above 
that 



C(n,c) = 0((Sc) 1/4 ) ■ exp 



where n,5 are as in (|46p . (|47[) . 



E 



•Xr, 



(52) 



In the following proposition, we describe another upper bound on |A n |, which is weaker than the 
one presented in Proposition [31 but has a simpler form. It is proved in Theorem [23] in Section [ 



Proposition 4. Suppose that n > is an even integer number, and that X n is the nth eigenvalue 
of the integral operator F c , defined via §5§ in Section POl Suppose also that 



2c /77T 

n > h V42. 

7T 



(53) 



Suppose furthermore that the real number rj(n, c) is defined via the formula 
T)(n,c) = 

V- (n + l)\ 7 



18 • (n + 1) 



• exp 



fa' F 



Xn - C z 



Xn 



Xn 



(54) 



where Xn is the nth eigenvalue of the differential operator L c , defined via ([3]) in Section^ and F, E 
are the complete elliptic integrals, defined, respectively, via (|39p . (1401) in Section \2.3[ Then, 



|A n | < rj(n,c). 



(55) 
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Remark 4. According to Proposition^ 



77(71, c) — 0(c) ■ exp 



as long as n is proportional to c. 



E 



(56) 



Both £(n, c) and n(n, c), defined, respectively, via (l50|) in Proposition[3]and (l54l) in Propositionl24l 
depend on Xni which somewhat obscures their behavior. In the following proposition, we eliminate 
this inconvenience by providing yet another upper bound on |A n |. The simplicity of this bound, as 
well as the fact that it depends only on n and c (and not on Xn), make Proposition [5] the principal 
result of this paper. 

It is proved in Theorem [32] in Section 14.31 and is illustrated via Experiment 3 in Section [5] 



Proposition 5. Suppose that c > is a real number, and that 

c > 22. 

Suppose also that S > is a real number, and that 



(57) 



3<*<-. 

Suppose, in addition,that n is a positive integer, and that 



2 

n > —c 



S ■ log 



' 4e7rc\ 



V s 



Suppose furthermore that the real number £(n, c) is defined via the formula 



£ (n, c) = 7056 • c • exp 



-SI 



s 

2nc 



The 



(58) 



(59) 



(60) 



|A„| < £(n,e). 



(61) 



3.2 Accuracy of Upper Bounds on |A n | 

In this subsection, we discuss the accuracy of the upper bounds on |A„|, presented in Propositions [3] 
SI [5] In this discussion, we use the analysis of Section |U previously reported results; and numerous 
numerical experiments, some of which are described in Section [5] Throughout this subsection, we 
suppose that n is a positive integer in the range 

— <n< — + 0(log(c)). (62) 

7T 7T 



According to the combination of Theorem [S] in Section 12.11 and Remark [21 

0(c 3/4 ), (63) 



|An| 
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where £(n, c) is that of Proposition^ On the other hand, both |A„ | and ((n, c) decay with n roughly 
exponentially at the same rate. Thus, the inequality (I5ip in Proposition [3] is reasonably tight (see 
also Experiment 2, Experiment 3 in Section [5]). 

The factor 0(c 3 / 4 ) in dMJ) is an artifact of the analysis in Section 14.11 The first source of 

inaccuracy is the inequality (|80p in the proof of Theorem 1111 In this inequality 

from above by 1, while numerical experiments indicate that 



(n,c) 



bounded 



< O^" 1 / 2 ), (64) 



for all integer k > 0. This contributes to the factor of order c 1 / 2 in (|6"3")l . The second source of 
inaccuracy is Theorem 1 141 which gives rise to the factor 

3 ' Xn/C 2 - 1 

in ([50)1 (see also Proposition [5]). This contributes to another factor of order c 1 / 4 in ([B"3"l) . 

In Propositions El [5] we introduce two additional upper bounds on |A„|, namely rj{n,c) and 
£(n, c). Due to Remarks [3j [4] and Proposition [5j 

77(n,c)=CKc)-0( C 3 / 4 ), 

£(n,c)=C(n,c)-0(c 3 / 4 ). (66) 

Thus, (|182j) is a tighter upper bound on |A n | than both (I19ip and (|226p . This is not surprising, since, 
due to Theorems 1241 [3"2l rj(n,c) and £(n,c) can be viewed as simplified and less accurate versions 
of £(n, c). There are two sources of the discrepancy (l66|) . First, in the proofs of Theorems l24l l32l 
the term (x„ — c 2 ) 1 ^ 4 is bounded from above by (^(c 1 / 2 ), while, in fact, it is of order c 1 / 4 (see (|6"5j) 
above). Additional factor of order c 1 / 2 in (|66p is due to Theorem |H] and Remark [2] in Section |2~T1 
See also results of numerical experiments, reported in Section [5] 

Finally, we observe that the upper bound v(n, c) on |A n |, introduced in Theorem2]in Section [2. 11 
is useless for n as in (|62p . due to the combination of Theorem IMl and Remark ITTI in Section |4~51 On 
the other hand, v(n, c) can be used to understand the behavior of |A„| as n — > oo, for a fixed c > 0. 



4 Analytical Apparatus 

The purpose of this section is to provide the analytical apparatus to be used in the rest of the paper. 
This principal results of this section are Theorems [231 [Ml 

4.1 Legendre Expansion 

In this subsection, we analyze the Legendre expansion of PSWFs, introduced in Section 12.21 This 
analysis will be subsequently used in Section [4.21 to prove the principal result of this paper. 

The following theorem is a direct consequence of the results outlined in Section 12.11 and Sec- 
tion [221 

Theorem 10. Suppose that c > is a real number, and n > is an even positive integer. Suppose 
also that the numbers a!^ ,c \ (4™' c \ . . . are defined via the formula 

4"' c) = J i> n (t) ■ P 2k -2{t) dt, (67) 
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for k — 1,2,..., where ip n is the nth PSWF corresponding to band limit c, and P k is the kth 
normalized Legendre polynomial. Then, the sequence |a[™' C ' ) | satisfies the recurrence relation 

a ■ 4"< c) + bx ■ a[ n ' c) = 0, 

Cfc+i • a k+2 J + b k+ i ■ a k+ { +c k - a\ = 0, (68) 
for k > 1, where the numbers c%, C2, . . . are defined via the formula 

2fc-(2fc-l) 2 

Cfc = ; V 1 • c 2 , 69 

(4fc-l) ■y/(Ak-3)-(Ak+l) 

for k > 1, and t/ie numbers b\, b%, . . . are defined via the formula 

» - 2 - » - 1) ■ (M - 1) + 2 • ^i').'^')- 1 ■ * - (TO 
/or fc > 1. i^ere % n i/ie n£/i eigenvalue of the prolate differential equation (|15p . Moreover, 

oo 

^w = E4 n,c) -^^(*) ! (71) 



,(«>c) 
fc=l 

and 

2 



E(4"' C) ) =1- (72) 



fe=l 



Proof. To establish (55} and t[7T]). we combine (|3"0]). ([55]). ([3"lf in Section O with Theorem [T] in 
Section 12.11 The identity (1721 follows from the fact that the normalized Legendre polynomials 
constitute an orthonormal basis for L 2 [— 1, 1]. ■ 

In the rest of the section, c > is a fixed real number, and n > is an even positive integer. 



The following theorem provides an upper bound on 
sequence. 



An,c) 



in terms of the elements of another 



Theorem 11. Suppose that the sequence ax, ot2, ■ ■ ■ is defined via the formula 



= "fa» (73) 



a 



(n,c) 
X 



for k > 1, where a^ 1 ' ^, a^"' "' , . . . are defined via (|67p in Theorem ] 1 (A Then, the sequence ax,ct2, ■ ■ ■ 
satisfies the recurrence relation 

ax = 1, 
a 2 = B , 

ctk+2 = B k ■ a k+1 - A k ■ a k , (74) 
for k > 1, where the sequence Ax, A%, . . . is defined via the formula 



fc-(2fc-l)-(4fc + 3) / 4fc + 5 
k [k + 1) • (2k + 1) • (4fe - 1) V4fc-3' 1 J 
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for k > 1, and the sequence Bo, B\, . . . is defined via the formula 



B k = 



Xn - 2k • (2k + 1)\ (4fc + 3) ■ V(4fc + 1) ■ (4k + 5) 

(2k + 1) • (2k + 2) 



(4k ■ (2k + 1) - 1) • y/(Ak + l) ■ (4k + 5) 
(4k -1) ■ (2k + 1) ■ (2k + 2) : 



for k > 0. Moreover, for every k — 1,2, 



< 



(76) 



(77) 



Proof. Due to (|68p in Theorem ITUl the recurrence relation ([74]) holds with Ak,Bk's defined via the 
formulae 



Cfe+l Ck+l 

where c k , b k s are defined, respectively, via (159"]) and (|7D|) . We observe that 



1 _ (4k + 3) • y/JW+T) ■ (4k + 5) 1 



Cfe+i (2fe + 1) • (2fc + 2) 

and readily obtain both ([75]) and ([76]). Next, due to ([72]) and ([73 

,(™, c ) 



1 > 



,(™. c ) 



,(™. c ) 



(n,c) 



(78) 



(79) 



(80) 



for all k = 1, 2, . . . , which implies ([77]). ■ 

It is somewhat easier to analyze a rescaled version of the sequence {a k } defined via (|73p in 
Theorem [IT] This observation is reflected in the following theorem. 

Theorem 12. Suppose that the sequence fi\,fii, • ■ • is defined via the formula 

Pk = a k 



4fc-3' 



(81) 



for k > 1, where a>i, a.2, . . . are defined via ([73]) m TTieoremQT] above. Suppose also that the sequence 
Bq, B\, ... is defined via the formula 



BX _ (4fc + l)-(4fc + 3) 



X„-c 2 -2fc- (2fc+l) 



fc (2fc+ 1) ■ (2fe + 2) 
/or fc > 0. Then, the sequence /?2, • ■ • satisfies the recurrence relation 



(82) 



h = V2, 

P 2 = B - s/2, 

Pk+2 — B k ■ /3fc+i — Ak ■ (3k, 



(83) 
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for k > 1, where Aq,A\,... are defined via the formula 

k- (2k- 1) • (4fc + 3) 



A, 



(k + 1) ■ (2k +1) • (4k- 1) 
/or fc > 0, and Bq, B\, . . . are defined via the formula 

B k = Bl + l + A k , 

for k > 0. 

Proof. Due to flU} and flU]), we have for all k = 1, 2, . . . 



(84) 



(85) 



2 


4fc + 


5 


4fc + 


1 


4fc + 


5 



• "fe+2 



4fc + 5 
i 



Sfc • a k+ i - 



4fc + 5 



and hence the recurrence relation 



4fc + 1 
holds with 



4fc- 3 
4fc + 5 



4fc- 3 



/4fc-3 , - /4fc + l 

Afc = ; - '4, B k = \ ——-B k . 



4fc + 5 



4fc + 5 



(86) 



(87) 



It remains to compute A k s and B k s. First, we observe that ([Mf follows immediately from the 
combination of (1751) with (l87l) . Second, we combine ((76)) with (|87|) to conclude that, for fc = 1, 2, . . . , 



"Xn - 2fc 


• (2k 


+ 1)1 


(4fc + 3) • (4fc + 1) 


(8k 2 4 


4fc- 


!)■ 


(4fc + 


1) 




c 2 




(2fe + 1) • (2fc + 2) (4fc-l) 


• (2fc 


+ 1 


)-(2fc 


+ 2) 


(Ak + 1) 


(Ak- 


f 3) 


~Xn - c 2 - 2fc • (2fc + 1)~ 


+ 










(2fc + l) 


(2k- 


f 2) 


c 2 












(4* + 3) 


■ (Ak 


+ !)• 


(4fc-l)-(4fc + l)-(8fc 2 + 4fc- 


-1) 












(4fc- 


1) • (2fc + l) • (2fc + 2) 












(Ak + 1) 


(Ak- 


f 3) 


~Xn - c 2 ~2k- (2k + 1)' 


+ 1 + 










(2fc + l) 


(2k- 


f2) 


c 2 










(Ak + 3) ■ 


(4fcH 


- 1) ■ (4fc - 1) - (4fc + l) • (8fc 2 


+ 4fc- 


1)- 


(4fc 


-1)- 


(2fc + l) 








(4fc - 1) • (2k + 1) • (2k 


+ 2) 








(Ak + 1) 


(4fc- 


f 3) 


'Xn-c 2 -2k-(2k+l)' 


+ 1 + 


A k , 








(2fc + 1) 


(2k- 


f 2) 


c 2 









which completes the proof. ■ 

The following theorem, in which we establish the monotonicity of both {a k } and {f3 k } up to a 
certain value of fc, is a consequence of Theorem 1121 

Theorem 13. Suppose that %„ > c 2 , and that /3i,/?2,--- are defined via (f5Tj) in Theorem \12\ 
Suppose also that the integer fcg is defined via the formula 



fc = max{fc = 1,2,... : 2fc ■ (2k + 1) < Xn - c 2 } 

ft 



= max < fc = 1, 2, . 



fc < 



Xn - C 2 + 



1 1 

4 ~ 4 



(89) 
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The 



V2 = Pi < 2 < ■ ■ ■ < Pko < Pk +i < Pk +2, 

and also, 

1 = oti < a>2 < ■ ■ ■ < a k() < a k()+ i < a ko+2 , 
where the sequences cmd {0k} ore defined via (|73[) and (|81[) . respectively. 

Proof. Due to (|85j) in Theorem [T2l and the assumption that x„ > c 2 , 

B = 1^ + 1>1. 



(90) 
(91) 



Therefore, due to (US in Theorem H2 



02 = B ■ fa > 



(92) 



(93) 



By induction, suppose that 1 < k < ko and assume that k < 0k+i- We observe that A k ,B k > 0, 
and combine this observation with (|53")) , (j54")l , (|85l) and ([59")) to conclude that 

Pk+2 = Pk+i + B k ■ p k +i + A k ■ (0 k+1 - k ) > j3 k+u (94) 

which implies (|90p . To establish (|9ip . we use (|81[) and observe that 



Oik 



4fc + l /? fe+1 /4fc + l 

> I t; r > 1, 



4fc - 3 k 



4fc-3 



(95) 



for all 1< k < k + 1. 



In the following theorem, we bound the sequence 0i,02, ■ ■ ■ , defined via (f5Tj) in Theorem [T2l by 
another sequence from below. 

Theorem 14. Suppose that \n > c 2 , and that the sequence pi,p 2 , ■ ■ ■ , is defined via the formula 



pk 



(4k - 6) • (4k - 4) • (Ak + 7) 
(4k - 2) • (Ak) ■ (Ak + 3) ' 

for all k — 1,2,.... Suppose also that the sequence A\ ew , A^" 1 , . . . 

Anew _ 4 _ 
^ifc — A k Pk, 

for all k — 1,2,..., where A k is defined via (|84[) in Theorem{T, 
sequence new , 0^ ew , ■ ■ ■ is defined via the formulae 



(96) 

is defined via the formula 

(97) 

J Suppose furthermore that the 



j^new 


= 01, 


onew 
P2 


= fa, 


onew 
P3 


= fo„ 


onew 
P/c+2 





for k > 2, where 0i,02 



are defined via (|81[). and Bi is defined via 



A new 
A k 



4k -4 4k - 6 4k + 7 
4k + 4 ' 4k + 2 ' 4k- 1' 



(98) 

in Theorem ] 12[ Then, 
(99) 
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for all k = 0,1, ... , and also 

= A n x ew < A% ew < A™ ew < ■ ■ ■ < Al ew < • • • < 1. (100) 

Moreover, 

/o anew . anew , , anew , anew , anew fir\l\ 

V2 = ft < ft <■■■ < P ko < Pk + 1 < Pfeo+2' ( 101 ) 

where fco is defined via (|89fl in Theorem ] 1S\ In addition, 

ft™<ft, ®T<h, ft 7i<^o+i, /3 fe "T 2 <Ao+2. (102) 



Proof. The identity (l99|) follows immediately from the combination of ((84)) and (f96|) . The mono- 
tonicity of {^4fe eu '} follows from the fact that, if we view Ak as a function of the real argument 

k, 

dA k _ (((3 + fc) ■ 8fc - 19) ■ 2fc - 51) ■ 8fc + 2 

~dfc~ ~ (4A:-l) 2 -(fc + l)2.(2fc + l) 2 ' ( ' 

which is positive for all fc > 2; combining this observation with the fact that A r k Lew tends to 1 as 
k —> co, we obtain (| 100[) . 

It follows from ([55]) by induction that ft^ > P"j+i as ^ on S as .B* > 0, which holds for all j < fc , 
due to (|S2J> and ([59]). This observation implies (fTOTj) . 

It remains to prove (|102[) . We observe that, due to (|96|) . the sequence = p\,p2,... grows 
monotonically and is bounded from above by 1. Combined with (|97[) . this implies that 

A n k ew <A k , k = l,2,.... (104) 

Eventually, we show by induction that 

fi+i ~ PT W < &+1 - Pk and ft^</3fc + i, (105) 

for all fc = 1, 2, . . . , fco + 1, with fco defined via (|B"9")) . For fc = 1,2, the inequalities (|105[) hold due to 
(JHHJ) - We assume that they hold for some fc < fc . First, we combine ([Ml), (ED), dFJJ), (EFJ, (|T(M|) and 
the induction hypothesis to conclude that 

- PT+i = B x k ■ PT+i + A r r • {PT+I - PT W ) < B x k • ft+i + A k ■ (ft +1 - ft). (106) 
Then, we combine ([52]). (|5T|) . (IMl) . dHFJ, (|104[) and the induction hypothesis to conclude that 

Pk+2 - m w 2 = ( B k + 1) • Wk+i - pt+i) + Ak • (ft+i - ft) - ai™ ■ - pt w ) 

> ft +1 - ft™ > 0, (107) 
which finishes the proof. ■ 

Theorem Q3] allows us to find a lower bound on ft by finding a lower bound on ft? e ™ , for all 
fc < fco + 2. In the following theorem, we simplify the recurrence relation (|98p by rescaling {(3 k ew }. 

Theorem 15. Suppose that Xn > c 2 + 6, and that the sequence ft 6 " 1 , ft ew , • ■ • is defined via (|M|) 
m Theorem \14\ Suppose also that the sequence /i , /a, . . . is defined via the formula 

(4fc-4)-(4fc-6) 
/fc 4fc^l ' (108) 
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for all k = 1, 2, . . . , and the sequence 71, 72, ■ • • is defined via the formulae 



anew 

71 = Pi > 

r anew 

Ik = Ik' P k , 

for k > 2. Then, the sequence 71,72, • ■ • satisfies, for k > 2, the recurrence relation 
7i = V2, 



72 



73 = 



7y/2 
16_yg 



• 2 + 3- 



Xn - C 



■ 3 + 15- 



Xn-C 2 105 Xn - C 2 Xn-C 2 - 6 105 

c 2 8 ' c 2 c 2 2c 2 



7fc+2 = (5fc+£f)-7fc+i-7fc, 
where the sequences {B^} and {-Bj^} are defined via the formulae 



Bi = 



4- (4fc + 1) • (4/c + 3) 2 
' k ~ 4k • (4k - 2) • (4k + 7) 

/or all k = 1, 2, ... , and 



Xn ~ c 2 - 2fc • (2k + 1) 



GO 



32fc 4 + 32fc 3 -38fc 2 + 7fc' 
for all k = 1,2, ... , respectively. Moreover, 

245 x« - c 2 - 6 _ , 



22 c 2 
where ko is defined via (|89p . and 

42 



B[>B I 2 >--->B{ > 0, 



^ = B* J > > • • • > Bi 1 > ■ ■ ■ > 2. 



(109) 

(110) 
(111) 

(112) 
(113) 

(114) 

(115) 

(116) 

(117) 



Proof. The identity (fTTTJ|) follows immediately from and (fTU9")) . Then, it follows from (|75 |) . (fTo) . 

that 



moreover, 



Bi = 



7V5 Xn-6 ll-y/5 7a/5 x« - c 2 - 6 7n/5 ll\/5 



12 4 

X« - c 2 - 6 



12 



We combine (|TT51) with ([73 ]) . (|8"I ]) . ([9"5 ]) . ([IDS]) . (fTTO)) to conclude that 

8 a 8 [2 8 [2 

72 = r /3 2 = --y-- Q2 = --y-- J B , 



(118) 



(119) 



(120) 
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from which (fTTTj) follows. Then we combine (fTTg)) . (fug]) with (0, (EH), ([Ml), (pM]) . (|TM|) to 

conclude that 



48 R 48 ^ 48 ^ rR , x 48x/2 

wv ~ 2 ( " ^ + 3.^vrio + 2i. ^-f- 6 vr 1; 



11 V 24 V c / V c / 6 

which simplifies to yield ([TT2"]) . The relation ([TT3"]) is established by using ([52]), (HHJ), (HZJ), (TTU5| . 
pUS)) to expand, for all fc > 2, 

7fe +2 = A+2 • = A+2 • {Bl + 1 + ^r 10 ) • PT+i - / fc +2 • Ar w ■ PT W 

= 7^ • {Bl + 1 + 4T) • 7 fc+i - ^ • AT W ■ l k . (122) 
Jfc+1 Jfc 

Since, due to dHZJ) , (HHFJ, we have 

fk+2 

Anew 

—? A k - 

Jk 

{An + 4) • (4n + 2) 4n - 1 {An - A) ■ {An - 6) • {An + 7) 

4n + 7 (4n - 4) • {An - 6) ' {An + 4) ■ (in + 2) • (4n - 1) ~~ ' ^ ' 

the identity (fTTB"]) readily follows from (fi"22"|) . (fT23)) . with 

^ = ^±2 . 73* (124) 



fk+ 

and 



B« = . {Al ew + l) . (125) 

Jfc+1 



We substitute (JB2>, (fTUS)) into ([T24]) to obtain (fTTlj) . Next, 



dA 7 



4- (4fc+ 1) ■ (4/c + 3) 2 



9 512 50 



4fc • {Ak- 2) • (4fc + 7) 
1 / 9 512 50 



14fc 2 21-(7 + 4fc) 2 3-(2fc-l) 



< 



(fc-1/2) 2 V 14 21-16 12y (A; — 1/2) 



< 0, (126) 



for all k > 1. Due to (l89l) . the term inside the square brackets of (|114[) is positive for all k > ko 
and monotonically decreases as k grows, which, combined with (|126[) . implies (|116p . Eventually, we 
substitute $7$, fTUS]) into (IT23|) and use (fl"2"3"]) to obtain, for all k > 1, 



ill fk+2 + Jk 



Jfc+1 



which yields (|115p through straightforward algebraic manipulations. The monotonicity relation (|117[) 
follows immediately from (|115p . ■ 



We analyze the sequence {jk} from Theorem [15] by considering the ratios of its consecutive 
elements. The latter are bounded from below by the largest eigenvalue of the characteristic equation 
of the recurrence relation (|113l) . In the following two theorems, we elaborate on these ideas. 
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Theorem 16. Suppose that \n > c > and that the sequence ri, ra, . . . is defined via the formula 

r k = —, (128) 

Ik 

for all k = 1,2, ... , where the sequence 71, 72, . . . is defined via (|109[) in Theorem ] 15[ Suppose also 
that the sequence o~i , 02, . . . is defined via the formula 



2 



.^ + fM, .,, (129) 



/or all k = 1,2,..., where B\, Bj. 1 are defined via (|114[) . (|115p in Theorem \lh\ respectively. Then, 

r 2 >S| + Bf. (130) 
Moreover, if B 2 + B 2 J > 2, iften ct 2 > 0, and 

r 2 > cr 2 . (131) 



Proo/. We use (flH| , (fTT5)) to obtain 



R ;, R H 44 121 X» - - 20 



(132) 



Next, we plug (|TTT|).([TT2l into (fT2"8")l to obtain 



7*2 = — • 3 + 15 



11 V c 2 8 c 2 c 2 2c 2 

2 + 3 '^^-) • ( 133 ) 
We subtract (|132|) from (|133p to obtain, by performing elementary algebraic manipulations, 

r (P. 1 1 - 247 1 1119 Xn ~° 2 98 (-2 1 o Xn-c 2 ^" 1 596 

247 98 398 

> 77-- 66 = 231 >0 ' (134) 

which implies (|130[) . Due to (|129D . a 2 is positive if and only if B 2 + B 2 ! > 2; in that case, 

B I 2 +B I 2 '>a 2 , (135) 

which, combined with (|130|) . implies (| 13 1[) . ■ 

The following theorem extends Theorem [TBI 

Theorem 17. Suppose that \n > c 2 , and that fco > 2, where fco is defined via (j89p in Theorem \13[ 
Suppose also that the sequences r% , r 2 , . . . and o~\, 02, • ■ • are defined, respectively, via (|128p , (|129[) 
in Theorem ] 16] Then, 

ai > a 2 > a 3 > ■ ■ ■ > a ko > 1. (136) 
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In addition, 

r 2 > r 3 > • • • > r ko > 1. (137) 

Moreover, 

r 2 > cr 2 > 1, r 3 > cr 3 > 1, . . . , r fco > cr fco > 1. (138) 



Proof. We combine ([Till) . (|TT51) . ([TTB]) . pT7|) in Theorem [T51 with (fHgf in Theorem \W\ to conclude 
that, for all k = 1, 2, . . . , ko, 

Bl + B k Bl 
o-fc > 2 > -y- > 1- (139) 

We use this in combination with (|116[) and (|117[) to conclude that (I136|) holds. Then, we use f|139[) 
and Theorem ITBl to conclude that 

r 2 > cr 2 > 1. (140) 

Next, we prove (|138l) by induction on k < ko. The case k = 2 is handled by (I140p . Suppose that 
2 < k < k , and (|138|) is true for k, i.e. 

r k > o-fe > 1. (141) 

We consider the quadratic equation 

x- 2 - (B, 7 + Bl 1 ) ■ x + 1 = 0, (142) 

in the unknown x. Due to (11 29[) and (I139[) . a k is the largest root of the quadratic equation (|142p . 
and, moreover, a^ 1 < 1 is its second (smallest) root. Thus, the left hand side of (|142[) is negative if 
and only if x £ (a^ 1 ,ak)- We combine this observation with (|141l) to conclude that 

rl^(B{ + B{ I )-r k + l>0, (143) 

and, consequently, 

r k> {B{ + B{ I )~-. (144) 

Then, we substitute (I128|) into (|113[) to obtain 

n+l = 2h±l = (Bl + Bl').>y k+1 ->y k = j + ^ J_ 

7fc+i 7fe+i " n- 

By combining fll44[) with (|145l) we conclude that 

r k > r k+1 . (146) 

Moreover, we combine (I141[) with (|145[) and use the fact that a k is a root of (|142[) to obtain the 
inequality 

r fc+1 = (i? 7 + O --> {B{ + ) - - = a k , (147) 
^fe o k 

However, combined with the already proved (|136[) and the fact that k < ko, the inequality (|147[) 
implies that 

r fc+ i>o- fc+ i. (148) 
This completes the proof of (I138p . The relation (I13T[) follows from the inequality (|146|) above. ■ 
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In the following theorem, we bound the product of several cr^'s by a definite integral. 

Theorem 18. Suppose that \n > c 2 , and that kg > 2, where ko is defined via (1891) in Theorem\ 
Suppose also that the real valued function g n is defined via the formula 



- 1 + 2 



Xn - <r 



2x 



\ 



1 + 2- 



Xn ~ C z 



2x 



(149) 



for the real values of x satisfying the inequality 4x 2 < \n — c 2 . Suppose furthermore that the sequence 
<7i,<72, • . . is defined via the formula (| 129[> in Theorem \16l Then, 



0-2 ■ 03 o-fco-i > (flVi(O)) 4 ' ex P 



(Vx^)/ 2 



log (#n 0*0) dx. 



(150) 



Proof. We observe that, for all k = 1, 2, . . . , 

4 • fc 2 < 2fc • (2ft + 1) < 4 • (fc + l) 2 < 2{k + 1) • (2(k + 1) + 1). 
In combination with (1891) . this implies that, for all k = 1, . . . , ko, 

Xn - c 2 - 4 • k 2 > 0. 

Moreover, due to (|114p . (|115l) in Theorem [T5l the inequality 



2 < 2 + 4 



Xn-c 2 (2-{k + \) 



< Bl + BU 



(151) 
(152) 

(153) 



holds for all k = 1, . . . , fco — lj where B^Bj, 1 are defined via (|114l) . (|115[) . respectively. We combine 
(1153[) with (|129l) in Theorem [16] and (|149l) above to obtain the inequality 

<T k >g n (k+l), (154) 

which holds for all fc = 1 , . . . , ko — 1 . Consequently, using the monotonicity of g n , 

o~2 ■ 0-3 Ofc _i > 

,o\ (AS (l S 9n{0) ■ ffn(l) 9n(k - 1) • ffn(fc ) 2 

9n 3 • .g„ 4 g n (k ) = -rr j- tt— > 

9n(0) ■ 3„(1) • g n {2) ■ g n (k ) 

g n (0y 4 ■ cxp (log(. 9n (0)) + • • • + log(. 9n (fc + 1)) + 2 • log(.g n (fc ))) 
Obviously, due to (|152p . the inequality 

rk+l 

log(g n (fc)) > J \og(g n (x)) dx 

holds for all k = 0, . . . , ko — 1. Next, due to (|8"9")l and (|151[) . we have 

1 



Therefore, 



&o < - Vxn - c 2 < k a + 2. 



2 ' log(5„(fc )) > -VXn - c 2 - ko ■ log(g n (k )) > / g„(x) dx. 



(155) 



(156) 



(157) 



(158) 



Thus, the inequality (|150p follows from the combination of (|155[) . (|156p and (|158p . 
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4.2 Principal Result 

In this subsection, we use the tools developed in Section 14.11 to derive an upper bound on |A n |. 
Theorem [22 is the principal result of this subsection. 

In the following theorem, we simplify the integral in (|150[) by expressing it in terms of elliptic 
functions. 

Theorem 19. Suppose that \n > c 2 , and that the real-valued function g n is defined via the formula 
HHg) in Theorem\M Then, 



log (aw (a;)) dx 



Xn - C 



2 r*/ 2 



sin 2 ((9) d6 



Moreover, 



L 



log(g n (x)) dx = y/x 



F 



E 



Xr, 



cos 2 ((9) 



>Xn_C_ 

Xn 



(159) 



(160) 



where F, E are the elliptic integrals defined, respectively, via the formula (|39[) , (140)) in Section 
Proof. We use Q149p and perform the change of variable 

2x 



VXr, 



(161) 



in the left-hand side of (|159l) to obtain 

•(Vx^)/2 



/ 

Jo 



log(g n (x)) dx 



\JXr, 



log g* 



ds 



V ■ 



■J log (l + 2F 2 (1 - s 2 ) + v/(l + 2V 2 (l - s 2 )) 2 - l) ds 



V-c 



log(h(s)) ds, 



where V is defined via the formula 



V = 



Xn - ^ 



and the function h : [0, 1] — > R is defined via the formula 



h(s) = 1 + 2^ 2 (1 - s 2 ) + y/(l + 2V 2 (1 -s 2 )) 2 -l. 
We observe that log(/i(l)) = and h(0) is finite, hence 



log(M«)) d» = [* • log(/i( a ))]J 



i ^s-h'is) 



h(s) 



ds = — 



1 s-h'(s) 
o ^( s ) 



ds. 



(162) 

(163) 
(164) 
(165) 



Then, we differentiate h(s), defined via (|164[) . with respect to s to obtain 

2- (1 + 2V 2 {1- s 2 )) ■ {-2V 2 -2s) 



h'(s) = -2V 2 ■ 2s + 
= -W 2 s • | 1 



2^(l + 2V r2 (l-s 2 )) 2 -1 

l + 2V 2 (l~s 2 ) \ 

y/{l + 2V 2 (l-s 2 )) 2 -l ) ~ v /(l + 21/ 2 (l- S 2 )) 2 -l' 



W 2 s ■ h{s) 



(166) 
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We substitute (|166[) into (|165p to obtain 



W 2 s 2 

los(hls)) ds = I — ; ds 
1 v/(l + 2V 2 (l-s 2 )) 2 -l 

AV 2 s 2 

ds 



o yJW A {\~ s 2 ) 2 +AV 2 {1- s 2 ) 

1 „2 



= 2V- = ds. (167) 

We perform the change of variable 

s = sin(0), ds = cos{9) ■ d9, (168) 

to transform (|167p into 

•io g( M.))^=^.r ,. si -. 2 . Wds (169) 



^l + l/ 2 - COS 2 (0) 

We combine (|162p . (|163D and (I169p to obtain the formula (|159[) . Next, we express f|159[) in terms of 
the elliptic integrals F(k) and E(k), defined, respectively, via (f39|) . (|40)) in Section l2~3l We note that 

r /2 k 2 sin 2 tdt k 2 r /2 sin 2 tdt , s 

F(k)-E(k)= = = - ■ . (170 



Motivated by (| 159[) and (|170[) . we solve the equation 

fc 2 _ Xn - c 2 
1 - k 2 ~ c 2 

in the unknown fc, to obtain the solution 



(171) 



k=J^—^-. (172) 



We plug ([T721 into (ITTO)) to conclude that 



F J^) - E J^) = ^ • r Sin2(g) ^ (173) 

VV xn J VV x» / <v%r y ^ - ^ . cos2 ^ 

We combine (fT5^|) with (|T73l to obtain ([Tip . ■ 

In the following theorem, we establish a relationship between the eigenvalue A n of the integral 
operator F c defined via Q in Section |2~T1 and the value of a± c ^ defined via (|67p in Theorem [TU] 

Theorem 20. Suppose that n > is an even integer number, and that X n is the nth eigenvalue of 
the integral operator F c defined via Q in Section \2.1l In other words, X n satisfies the identity §5§ 



in 



in Section \2.1[ Suppose also, that the sequence a\ ' , a\ ' , ■ ■ ■ is defined via the formula (I67[) 
Theorem \1(A Then, 

where i\) n is the nth prolate spheroidal wave function defined in Section \2.1\ 
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Proof. Due to © in Section O {25]), (JUJ in Section O and ([57]) above, 

X n -M0)= i>„(t)dt = V2- ip n (t)-Po(t)dt = V2-ap c) , (175) 



from which (|174l) readily follows. ■ 

In the following theorem, we provide an upper bound on |A n | in terms of the elements of the 
sequence {7fc}, defined via (|109l) in Theorem IT51 above. 

Theorem 21. Suppose that n > is an even integer number, and that X n is the nth eigenvalue 
of the integral operator F c , defined via ([3 ]) . ([5 ]) in Section \ 2.1\ Suppose also that \n > c 2 , &nd 
that ko > 2, where ko is defined via (|89p in Theorem ] 1S\ Suppose furthermore, that the sequence 
71,72,... is defined via (|109p in Theorem \l,5l Then, 

\\ I c 2 (4fco ~ 4) ■ (4fc - 6) \_ 

' U{ \M0)\ (4fc - 1) • V4fc^3 lk0 - 1 ' 

Proof. We combine the inequality ([77]) in Theorem [11] with the identity (|174[) in Theorem [551 to 
conclude that 

\\ I = — I ( ™' c) l <r — — 2 1 1 /, 77 x 

1 1^(0)1'' 1 1 < |^ n (0)| ' a feo |^(0)| ' V4fc^3'/?fco' 1 J 

where /3 fco is defined via ([SI]) in Thcorem[TJ Next, we combine 0, (1T52"]) in Theorem[U ([155]) , CE55]) 
in Theorem [15] and (|177p to obtain the inequality 

, 2 112 1 1 

|A„| < — -— • • — < 



|^n(0)| ^41^3 /3 feo " |V«(0)| V4fc^3 /3 fc " et " 
2 (4fc - 4) ■ (4fc - 6) 1 



(178) 



IV™ (0)| (4fc - 1) • V4fc - 3 7fc ' 
which is precisely (|176l) . ■ 

The following theorem is a direct consequence of Theorems [5J [Jj in Section 12.11 
Theorem 22. Suppose that n > is a positive integer. Suppose also that n > (2c/n) + \/42- Then, 

Xn>c 2 + 42, (179) 

and also, 

k > 2, (180) 

where ko is defined via (|89[) in Theorem \13\ 

Proof. Suppose that c 2 < \n > c 2 + 2. Then, due to Theorem [5J 



n< 2 - [\F^dt< 2 - [\l C i + ^dt 



v 1-t 2 -ttJo V 

2c 2Vl2 f 1 dt 2c r— , . 

<— + • / - — + V42. 181 

We combine (|181l) with Theorem [5J to conclude (|179[) . Then, we combine (|179p with ([55]) in Theo- 
rem [TJ to conclude (pT551) . ■ 
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The following theorem is the principal result of this paper. 

Theorem 23. Suppose that n > is an even integer number, and that X n is the nth eigenvalue of 
the integral operator F c , defined via (HJ), in Section [KT[ Suppose also that Xn > <? + 42. Suppose 
furthermore that the real number £(n, c) is defined via the formula 

7 ' ^ 



C) = oL; , ^ • , - — 2 — • (Xn - C ) 



2|V„(0)| 3- Xn /c 2 -l 



exp 



Xn - ^ 



- E 



Xr. 



X;, 



(182) 



where F, E are the complete elliptic integrals, defined, respectively, via (|39p , f)40[) m Section 
Then, 



|A„| < C(«,c). 



(183) 



Proof. We start with observing that, due to (181)1) in Theorem [TO] and (|157[) in Theorem [TO] the 
inequality Xn > c 2 + 42 implies that fc > 2. We combine ([1051) in Theorem ITOl (fHgj) . (fT2"0|) in 
Theorem [TO] and (|138[) in Theorem [T7] to obtain the inequality 



lk =12 



73 



7Ao-i 7fe 



72-^2 r feo _ 2 ■ r fco _x > 72 ■ (er 2 



72 7fc -2 7*o-i 

Next, we substitute (|149l) . (|150l) in Theorem [TBI into (|184l) to obtain the inequality 

-(\/x^)/ 2 



°fco-iJ 



(184) 



7fc > 72 ■ (Sn(0)) • exp 

Xn - I 



> 72 • 2 + 4 



exp 



log(g„(x)) 



log(£f„(a;)) cte, 



(185) 



where the function g n is defined via (|149p . Then, we plug the identity (|159[) from Theorem [TO] into 
(|185p to obtain the inequality 



1 1 

7*o 72 



2 + 4 



Xr 



■ exp 



c 2 f n/2 sm 2 (0)d0 



cos 2 ((9) 



We use (J50J in Theorem H51 and (fTOTj) in Theorem H51 to conclude that 



(4fc - 4) • (4fc - 6) 



.2\ 4 



(4fc - 1) • V4fc - 3 
We substitute <[T57|) into (JTT0|) in Theorem [211 to obtain 



IAJ < 



V2- 



[Xn - C 



,2^ i 



1 



' ' l^n(0)| ' ' lk 

Finally, we combine (jllll) in Theorem H5l with (|186l) . (11881) to obtain 

JN "I 



(186) 



(187) 



(188) 



|A„| < 



7 



21^(0)1 



(Xr, 



„2\ 4 



2 + 3- 



Xn - cr 



2 + 4- 



Xn - cr 



exp 



n/2 



sin 2 (6») d0 



cos 2 ((9) 



(189) 
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Eventually, we combine (I160p in Theorem [H] with fl!89[) to conclude (|183p . ■ 
Remark 5. The assumptions of Theorem \23\ are satisfied if n is an even integer such that 

n>— + V42, (190) 

since, in this case, Xn > <? + 42 due to Theorem\22\ 

4.3 Weaker But Simpler Bounds 

In this subsection, we use Theorem |2"31 in Section 14^21 to derive several upper bounds on |A„|. While 
these bounds are weaker than £(n, c) defined via (|182[) . they have a simpler form, and contribute to a 
better understanding of the decay of |A„|. The principal results of this subsection are Theorems [Ml 

M 

In the following theorem, we simplify the inequality ([1831) . The resulting upper bound on |A n | is 
weaker than (|183[) in Theorem [531 but has a simpler form. 

Theorem 24. Suppose that n > is an even integer number, and that A„ is the nth eigenvalue of 
the integral operator F c , defined via ([5J) in Section \2.1\ Suppose also that \n > c 2 +42. Suppose 
furthermore that the real number 77(71, c) is defined via the formula 



n(n,c) = 18 • (n + 1) 



7r- (n + l)\ 7 



exp 



X^- F J^— -E 



Xn \ V Xn 



(191) 



where F, E are the complete elliptic integrals, defined, respectively, via (I3U1) . fl4"D]) in Section \2. 3[ 
Then, 

\X n \<n(n,c). (192) 



Proof. We use ((23)) in Theorem \8\ in Section |2~T1 to conclude that 

(Xn-c 2 ) 1/4 <(x«) 1/4 <(f -(n + l)) 1/2 . (193) 

Next, 

(2 + 3 • ^) _1 • (2 + 4 • ^) 4 < 2? • ) 3 . (194) 
We combine Theorems E] in Section O with ([T9"5)i . (fTM)) to conclude that 
1 2i< 



\M0)\ 3- X „/c 2 -l 

(4- X »/c 2 -2) 

3-Xn/c 2 -l 



4. (B + 1) V..jr.(^)'".(E. (B + 1) ) 



7/2 /7r _ N l/2 

< 



/tt 7 / % /rr-fn+nA 7 /F , . /7r-(n+l)\ 7 , , 

We conclude by combining the inequality (|183[) in Theorem above with the inequality (|195l) . ■ 
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Both C(n, c) and 77(71,0), defined, respectively, via (|182l) in Theorem [231 and (|19ip in Theorem |2"H 
contain an exponential term (of the form exp [...]). This term depends on band limit c and prolate 
index n through \ n , which somewhat obscures its behavior. The following theorem eliminates this 
inconvenience. 

Theorem 25. Suppose that n is a positive integer such that n > Icjix , and that the function 
f : [0, 00) — > M is defined via the formula 

rr/2 



f{x) = -l + y/x + cos 2 (6) d0. (196) 



Suppose also that the function H : [0, 00) — > M is the inverse of f , in other words, 

y = f(H(y)) = -l+ y/ H(y)+ cos* (0)M, (197) 
Jo 

for all y > 0. Suppose furthermore that the function G : [0, 00) — > K is defined via the formula 

r* (198) 

Jo yl + x ■ cos z [p) 



for all x > 0. Then, 



Moreover, 



«■ * (i - • - (* (i - o) < ^- (* (f^j - * (v^j) • « 

where F,E are the complete elliptic integrals, defined, respectively, via (|39[) . (|40p m Section \2.tA 

Proof. Obviously, the function /, defined via (|196[) . is monotonically increasing. Moreover, /(0) = 0, 
and 

lim f{x) = 00. (201) 

x — >oo 

Therefore, i?(j/) is well defined for all y > 0, and, moreover, the function is monotonically 
increasing. This observation, combined with Theorems [5J E] in Section 12.11 implies the inequality 

(HMD. 

Next, the right hand side of (|200j) increases with \n, due to the combination of (j3"9")l . (|40l) in 
Section 12.31 This observation, combined with (|173[) in the proof of Theorem \T9\ (|198p and (|199p , 
implies POUI) . ■ 

Remark 6. TTie functions H,G, defined, respectively, via (|197|) . (|198p above, do not depend on 
either of n,c,\ n . Therefore, while the right-hand side of (|200p does depend on %n> ^ s left-hand side 
depends solely on c and n. 

In the following theorem, we provide simple lower and upper bounds on H, defined via (|197p in 
Theorem 
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Theorem 26. Suppose that the function H : [0, oo) — > R is defined via (|197p in Theorem \25\ Then, 

s < H (i.\og l ^f\< S + S j, (202) 

for all real < s < 5. 

Proof. The proof of (|202j) is straightforward, elementary, and is based on Pdj) in Section 12.31 it will 
be omitted. The correctness of Theorem [55] has been validated numerically. ■ 



Remark 7. The relative error of the lower bound in (|202j) is below 0.07 for allO < s < 5; moreover, 
this error grows roughly linearly with s to w 0.0085 for all < s < 0.1. The relative error of the 
upper bound in (|202l) grows roughly linearly with s to 1, for all < s < 5. 



In the following theorem, we provide simple lower and upper bound on G, defined via (|198p in 
Theorem [25] 

Theorem 27. Suppose that the function G : [0,oo) — > R is defined via (|198p in Theorem \25\ Then, 

J. (l- |) (203) 

for all real < x < 5. 

Proof. The proof of (|203[) is elementary, and is based on the fact that, for all x > 0, 

G(*) = J-(l-f + ^ + 0(*»)), (204) 

where G is defined via (1198[) in Theorem [25] The correctness of Theorem [27] has been validated 
numerically. ■ 

Remark 8. The relative errors of both lower and upper bounds in (|203p are below 0.6 for all 
< x < 5; moreover, these errors are below 0.01 for all < x < 0.1, and grow roughly linearly with 
x in this interval. 

The following theorem is in the spirit of Theorems [2T1 

Theorem 28. Suppose that the functions H,G : [0,oo) — > R are defined, respectively, via (|197p . 
(TTM1) in Theorem\M Then, 

/or a/Z reaZ < s < 5. Moreover, the function x — > H(x) ■ G{H(x)) is monotonically increasing. 

Proof. The proof is based on Theorems 1261 1271 is elementary, and will be omitted. The correctness 
of Theorem |2"81 has been validated numerically. ■ 

Remark 9. The relative errors of both lower and upper bounds in (|205p are below 0.5 for all 
< s < 5. Moreover, these errors are below 0.01 for all < s < 0.1, and grow roughly linearly with 
s in this interval. 

The following theorem is a consequence of Theorems [55] - 
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Theorem 29. Suppose that 8 > is a real number, such that 



4 

Suppose also that n is a positive integer, such that 



0<6<^--c. (206) 



n>-c+\-5-\og( A -^). (207; 

7T TT Z \ 



Then, 



where F,E are the complete elliptic integrals, defined, respectively, via (|39l) . (|40[) m Section \2.tA 
Proof. It follows from ((207)) that 

^^^•c' 108 ^)- (209) 
We define the real number s > via the formula 

s = — , 210 

7TC 

and observe that < s < 5 due to POSI) . We combine (12"W]) . (I2TU)) and Theorem [25] to obtain 

«(i-')- G (-(i- 1 ))> 

H (iAi og (^y G ( fl (i.:u g ^ 



7T C \ S J J \ \7T C \ (5 

We substitute (|2TTj) into the inequality ([2"0U1) in Theorem [25] to obtain (|2H5|) . ■ 

In the following theorem, we derive an upper bound on \n in terms of n. 
Theorem 30. Suppose that n is a positive integer, and that 



for some 



4 

Then, 

Xn-c 2 8 S 



2c 2c 2 p , /4evrc\ „ 

- < n < - + -j • 5 • log — — - 3, (212) 



3 < 5 < — ■ c. (213) 



< . (214) 

7T C 
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Proof. We combine QUI]) , (J5T3[) . (ITM| in Theorem [55] and ([50"5]l in Theorem [25] to obtain 

Xn ~ C 2 / (n + 3) ■ 7T 













-a 


/ 7 < 7TC 







. V <lf U'c- l0S l- 

which implies ([214[1 . We also observe that (|2 1 3[) implies 

2 r , /4e?rc\ 
- 5 .J.log(- r -j-8>1.3 > 

and hence there exist integer n that satisfy (|212p . 

In the following theorem, we derive an upper bound on the non-exponential term of £( 
defined via ([Tg5jl in Theorem [53] 

Theorem 31. Suppose that n is an even positive integer, and that 

2c 2 2 . . /ienc\ 

— < n < -c + — ■ S ■ log ( — ^ j - 3, 



V S 



for some 



Then, 



3 < 5 < — • C. 

4 



2|V„(0)| 3-Xn/c 2 
6 \n J \7rc/ \ 7rc 



Proof. We use ([214]) to obtain 



(4- X „/c 2 -2) 4 _4 (4.( X „-c 2 )/c 2 + 2) £ - 



< 32Y 1 + 2 .X ! ^y < 32 ^ , 16^ 3 



3 • xn/c 2 - 1 3 4- ( Xn - c 2 )/c 2 + 8/3 

" 3 

Then, we use (|214[) to obtain 

<--^<(^) > - 

Next, we combine Theorems [7] [9] in Section |2~T1 with Theorem l30l to obtain 

|Vn(0)| V c C C \ TTC J \ 7TC 

We combine Q550]) , (E2U), ([5521 to obtain ([51511 . 

The following theorem is one of the principal results of this subsection. 
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Theorem 32. Suppose that c > is a real number, and that 

c > 22. 

Suppose also that 5 > is a real number, and that 

3<5< ~6- 

Suppose, in addition, that n is a positive integer, and that 

2c 2 „ /4e?rc\ 

Suppose furthermore that the real number £(n, c) is defined via the formula 

5 



£(n, c) = 7056 • c • exp 



5 1- 



2nc 



Then, 



| A„ | < £(n,c). 



Proof. Suppose first that n is an even positive integer of the form 

2c 2 r , /4e7rc 
n = — + — • <5 • log — — 
ir ir z \ 6 



for some 3 < 5 < 7rc/16 (in other words, (|225p is an identity rather than an inequality) 
that, for all real t > 0, 



dt 



t ■ log 



(nr)) 



log 



/ 47TC 



We combine ([2231) with (I229| to obtain 

4e7rc 



2 /7TC 

-5 • "5- ' lo § 



7T^ 



M/8 



5 -lot 



/ 4e7rc\ 



2 / 7TC A / 4-7TC 

> ■ (t - *) ' log l 

> ^- log (32) >3. 

87T 



Therefore, it is possible to choose a real number 5 such that 

3 < 5 < — , 



and also 



2c 2 » /4 e7 rc\ 
ri=- + — log f-p J -3. 



7T 7T 

Due to the combination of (1231[) , (12321) and Theorem [2U 

7 (4- Xn /c 2 -2^ 



2|^(0)| 
448 



• Xn/C 2 - 1 



32 

1+lgl = 7,):,(i 
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We observe that the right-hand side of (|233[) is independent of S. We combine this observation with 
(12331) , (|183[) in Theorem [231 (|208p in Theorem [311 and the fact that | A n | decrease monotonically 
with n, to obtain (|227[) . ■ 



Definition 1 (S(n)). Suppose that n is a positive integer, and that 



2c 10c 
— < n < . 

7T 7T 



We define the real number 5(n) to be the solution of the equation 



2c 2 
n = 1 - • 

ir 7r z 

the unknown X in the interval < X < 4nc. 



X-log 



Aenc 
X 



(234) 



(235) 



Remark 10. We observe that the right-hand side of (I235P is an increasing function of X in the 
range < X < Awe, due to (|229[) in the proof of Theorem \33{ Therefore, S(n) is well defined. 

In the following theorem, we derive yet another upper bound on |A n |. 

Theorem 33. Suppose that n > is a positive integer, and that n > (2c/ir) 
that the real number defined via the formula 

Xn 

•En — n~ • 
C z 

Then, 



42. Suppose also 
(236) 



|A„| < 1195 -c- (x n )* ■ (x n - 1)3 



• exp 



7T 

1 



(237) 



Proof. We use l|236|) to obtain 

(4. X »/c 2 -2) 4 _4 (4 ■ (xn c 2 )/c 2 + 2) 4 256 / _1\ 3 
3-Xn/c 2 -l 3 ' 4 • ( X „ - + 8/3 K 3 ' [ Xn 2j ' [Z ™> 

Next, we combine Theorems El [§] in Section |2~T1 and (|236[) to obtain 

' <4-v^- A /^< --{XnY* = 4>/c ■(*»)*■ (239) 



|Vn(0)| V c 2 c 

We combine flUS} and {255)1 to obtain 



7 (x„~c 2 )* (4- X „/c 2 -2) 3 i / 1\ 

<1195-c-(x n )*-(af„-l)*- hc n -- . (240) 



2 |^ n (0)| 3- Xn /c 2 -l 
Also, we combine fl35J|, (gDJl in Section [221 with (j!70|) in the proof of Theorem [Hfl to obtain 



' Xn - C 2 \ / /Xn - C 2 \ 7T Xn 



F fir- S V ri'V (241) 



We combine (|356l . (jMJl) . (|MT1) with Theorems |2"21 |2"31 to obtain (g57)) . 
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We conclude this subsection with the following theorem, that describes the behavior of the upper 
bound v(n,c) on |A„| (see (JT7J) in TheoremHin Section 133]) . 



Theorem 34. Suppose that n is a positive integer, and that 

1 

25 



- • c < n < I - + — ) • c. (242) 



The 



v(n,c)>± t (243) 
where v(n,c) is defined via (|16[) in Theorem^in Section \2.1{ 

Proof. We carry out elementary calculations, involving the well known Stirling's approximation 
formula for the gamma function, to obtain the inequality 

\f2~wn / ce \ ™ . 

vin, c) > — , (244 

V ' J ~ 2n+l V4n/ ' V ' 

for all n in the range (I242p . We use (|244j) to obtain the inequality 

1 , / ce ' 



1 / ce \ 

log(i/(n,c)) >log— =+n-log(^— J 



>---log(c)+ - + — -c-log' 



2 ow V 71 " 25 / V 2 / 7r + !/ 25 

> -i • log(c) + > i . (1 - log(250)) > -2.27. (245) 



The inequality ((243)) follows directly from (|245j) . ■ 

Remark 11. According to Theorem \34\ the inequality (|17[) o/ Theorem^ in Section \2.1\ is trivial 
for all integer n < (2/ir + 1/25) • c. In particular, for such n this inequality is useless. 



5 Numerical Results 

In this section, we illustrate the results of Section U via several numerical experiments. All the 
calculations were implemented in FORTRAN (the Lahey 95 LINUX version) and were carried out 
in either double or quadruple precision. The algorithms for the evaluation of PSWFs and the 
associated eigenvalues were based on [3]. 



5.1 Experiment 1 

In this experiment, we demonstrate the behavior of |A n | with < n < 2c/~k, for several values of 
band limit c > 0. 

For each of five different values of c = 10, 10 2 , 10 3 , 10 4 , 10 5 , we do the following. First, we 
evaluate |A„| numerically, for n = 0, n « c/-k and n « 2c/n. For each such n, we also compute 

= ( c / 27r ) ' |A n |. Here A n is the nth eigenvalue of the integral operator F c , and /i n is the nth 
eigenvalue of the integral operator Q c (see (jU), ([5]), ©, (flUl) in Section l2~Tj) . 

In addition, we fix c = 100, and evaluate |A n | numerically, for all integer n between and 2c/ir. 

The results of Experiment 1 are shown in Table []] and Figure [T] Table Q] has the following 
structure. The first two columns contain the band limit c and the prolate index n, respectively. The 
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c 


n 


M/(2c) 


|Anl 

1 " 1 


/i„ = (c/2tt) 


|A„| 2 


10 





0.00000E- 


-00 


0.79267E+00 


0.10000E+01 


10 


3 


0.47124E- 


-00 


0.79183E+00 


0.99790E- 


-00 


10 


6 


0.94248E- 


-00 


0.52588E+00 


0.44015E- 


-00 


100 





O.OOOOOE- 


-00 


0.25066E+00 


0.10000E- 


-01 


100 


31 


0.48695E- 


-00 


0.25066E+00 


0.10000E- 


-01 


100 


63 


0.98960E- 


-00 


0.18589E+00 


0.54997E- 


-00 


1000 





O.OOOOOE- 


-00 


0.79267E-01 


0.10000E- 


-01 


1000 


318 


0.49951E+00 


0.79267E-01 


0.10000E- 


-01 


1000 


636 


0.99903E- 


-00 


0.57640E-01 


0.52877E- 


-00 


10000 





O.OOOOOE- 


-00 


0.25066E-01 


0.10000E- 


-01 


10000 


3183 


0.49998E- 


-00 


0.25066E-01 


0.10000E- 


-01 


10000 


6366 


0.99997E- 


-00 


0.16644E-01 


0.44088E- 


-00 


100000 





O.OOOOOE- 


-00 


0.79267E-02 


0.10000E- 


-01 


100000 


31830 


0.49998E- 


-00 


0.79267E-02 


0.10000E+01 


100000 


63661 


0.99998E- 


-00 


0.60295E-02 


0.57861E- 


-00 



Table 1: Behavior of |A„| for < n < 2c/ tt. Corresponds to Experiment 1 in Section^ 

third column contains the ratio of n to 2c/tt. The fourth column contains |A„|. The last column 
contains the eigenvalue fi n of the integral operator Q c (see ([6]), (|10[) in Section l2TTj) . 

In Figure [I] we plot |A n |, corresponding to c = 100, as a function of n, for integer n between 
and 2c/n. 

Several observations can be made from Table [T] and Figure [TJ 

1. For all five values of band limit c, the eigenvalue fi n decreases from « 1 to « 1/2, as n 
increases from to (2c/ w). In other words, the first 2c/ V eigenvalues A„ have roughly the 
same magnitude w ^j2%/c. This observation confirms Theorem [5] in Section 12.11 

2. Due to Theorem|6]in Section [2~T1 the bounds on the decay of |A„|, established in Section^ hold 
for n greater than 2c/7r only (see also Remark[S]). Thus, Tabled] indicates that this assumption 
on n is, in fact, not restrictive, since the first 2c/ tt eigenvalues have roughly constant magnitude. 

5.2 Experiment 2 

In this experiment, we illustrate Theorem 1231 As opposed to Experiment 1, we demonstrate the 
behavior of |A n | for n > 2c/n. 

In this experiment, we proceed as follows. First, we pick band limit c > (more or less arbitrar- 
ily). Then, for each even integer n in the range 

— < n < — + 20 ■ log(c), (246) 

7T TT 

we evaluate numerically |A„| and £(n, c), where the latter is defined via (|182[) in Theorem [231 

The results of Experiment 2 are shown in Figures [5] - @] and in Tabled] In Figures [5] - 2] we plot 

both log(|A„|) and \og(((n, c)) as functions of n. Each of Figures [5] -0] corresponds to a certain value 

of band limit (c = 10, 10 2 , 10 3 , 10 4 , 10 5 , respectively). 

Table [5] has the following structure. The first column contains precision e = e~ 50 ,e~ 100 . The 

second column contains band limit c. The third column contains the integer nx(e), defined via the 
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ni(e) 


Ai(e) 


"2(e) 


A 2 (e) 


n 2 {e) - m(e) 


e 


-50 


10 


32 


0.11133E- 


-02 


38 


0.13738E+02 


6 


e 


-50 


10 2 


107 


0.94107E+01 


114 


0.10931E+02 


7 


e 


-50 


10 3 


700 


0.91752E- 


-01 


712 


0.10912E+02 


12 


e 


-50 


10 4 


6450 


0.90987E- 


-01 


6468 


0.11053E+02 


18 


e 


-50 


10 5 


63765 


0.89484E- 


-01 


63792 


0.11294E+02 


27 


e~ 


-100 


10 


50 


0.18950E- 


-02 


56 


0.21556E+02 


6 


e" 


-100 


10 2 


138 


0.16142E- 


-02 


146 


0.17879E+02 


8 


e" 


-100 


10 3 


753 


0.16848E- 


-02 


764 


0.18440E+02 


11 


e" 


-100 


10 4 


6526 


0.17350E- 


-02 


6542 


0.19087E+02 


16 


e" 


-100 


10 5 


63864 


0.17547E- 


-02 


63890 


0.19806E+02 
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Table 2: Illustration of Theorem\23[ Corresponds to Experiment 2 in Section^ 



formula 



ni(e) — rnin {fc > 2c/tt : |Afc|<e}. 



In other words, ni(e) is the integer satisfying the inequality 

|A ni ( £ )_i| > e > |A ni ( £ )|. 



(247) 



(248) 



The fourth column contains Ai(e), defined to be the difference between 711(e) and 2c/ it, scaled by 
log(c). In other words, 



Ai(e) 



ni(e) — 2c/ir 



log(c) 

The fifth column contains the even integer 712(e), defined via the formula 
712(e) = min{fc > 2c/7r : fc is even, |£(fc, c)| < e} . 

k 

In other words, 77.2(e) is the even integer satisfying the inequality 

|C(n 2 (e)-2,c)| >e> \((n 2 (e),c)\. 



(249) 



(250) 



(251) 



The sixth column contains A2(e), defined to be the difference between 712(e) and 2c/7r, scaled by 
log(c). In other words, 



A 2 (e) 



77 2 (e) - 2c/7T 

log(c) 



(252) 



The last column contains the difference between 772(e) and 77i(e). 
Several observations can be made from Figures [5] - H] and Table [5] 

1. In all figures, |A n | < C( n , c), as expected, which confirms Theorem [231 

2. For each c, both |A n | and £(77, c) decay roughly exponentially fast with 71. 

3. For each c, both |A ra | and C( n J c ) decrease to roughly e~ 125 , as n increases from 2c/7r to 
2c/7r + 20 • log(c). In particular, 



A 



2c/7T+20-log(c) I 



,-125 



(253) 
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for c = 10, 10 2 , 10 3 , 10 4 , 10 5 . The fact that the right-hand side of ([2331 is the same for all c is 
somewhat surprising. However, this is not coincidental, as will be illustrated in Experiment 3 
below. 

4. For c = 10 2 , 10 3 , 10 4 , 10 5 , it suffices to take n w 2c/?r + 9 • log(c) to ensure that |A„| « e" 50 
(see third column in Tabled). In addition, it suffices to take n w 2c/ it + 17 ■ log(c) to ensure 
that | A« | ~ e~ 100 . In other words, 

Tii (e) « ^ + 0.17 • log • log(c), (254) 

where ni(e) is defined via (|247p above (see also (12531) ). 

5. For c = 10 2 , 10 3 , 10 4 , 10 5 , it suffices to take n 2c/?r + 11 • log(c) to ensure that C(«, c) « e~ 50 
(see fifth column in Tabled]). In addition, it suffices to take n ps 2c/7r + 19 • log(c) to ensure 
that C(n,c) « e~ 100 . In other words, 

„ a ( e )w ^ +O .2-log0J .log(c), (255) 

where 712(e) is defined via (|250[) above (see also (|253D . (|254jl ). 

6. The difference 712(e) — Tii(e) is roughly independent of e, and grows only slowly as c increases 
(see last column of Table [2]) • In other words, suppose that one needs to determine 71 such that 
|Afe| < e -50 for all k > n. Due to (|247[) . iii(e -50 ) would be the minimal such n. On the other 
hand, 71 = 7i2(e -50 ) is only larger by 6 for c = 10 and by 27 for c = 10 5 . 

5.3 Experiment 3 

In this experiment, we illustrate Theorem 1321 We proceed as follows. First, we pick band limit c > 
(more or less arbitrarily). Then, we define the positive integer 7i max to be the minimal even integer 
such that 

> - + \ ■ 150 • log (^) « - + 30.4 • log(0.23 • c). (256) 

TT 1T Z \ 150 / 7T 



'■max 



Then, for each positive even integer n in the range 



2c . s 

— < n < 7i max , (257) 

7T 



we evaluate the following quantities: 

• the eigenvalue A„ of the operator F c (see (J4j), (J5]) in Section l2~Tj) ; 

• <5(ti) of Definition Q] in Section |4~51 

• ({n,c), defined via ([T52l in Theorem [2^1 in Section fQ 

• £(n,c), defined via (|226D in Theorem I3"2l in Section |4~51 



The results of Experiment 3 are shown in Figures 5(a)| |5(b) that correspond, respectively, to band 
limit c = 10 4 and c = 10 5 . In each of Figures 5(a) 5(b) we plot log(|A„|), —S(n), log(C(n, c)) and 



log (£(n,c)) as functions of n. 

Several observations can be made from Figures 5(a)| 5(b) [ and from more detailed experiments 
by the author. 
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1. In both figures, 



log(|A„|) < -S(n) < log(C(n,c)) < log(£(n,c)), 



(258) 



for all n. This observation confirms both Theorem [531 of Section 14.21 and Theorem (33] of 
Section [4.31 Also, £(n,c) is weaker than £(n,c) as an upper bound on |A„|, as expected. 

2. All the four functions, plotted in Figures 5(a)[ 5(b) decay roughly exponentially with n. 
Moreover, 



log(|A„|) « log 



5{n), 



(259) 



in correspondence with Theorem[5]in Section [2.1l In particular, even the weakest bound £(n, c) 
correctly captures the exponential decay of |A n |. On the other hand, £(n, c) overestimates |A„| 
by a roughly constant factor of order c 3 / 2 (see also Section [ 
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Figure 1: Behavior of |A ra | for < n < 2c/ir, with c = 100. Corresponds to Experiment 1 in 
Section [31 
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Figure 2: Illustration of Theorem |2"31 with c = 10. Corresponds to Experiment 2 in Section [5j 
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(b) c = 1,000. 



Figure 3: Illustration of Theorem |2"51 Corresponds to Experiment 2 in Section [5] 
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Figure 4: Illustration of Theorem |2"51 Corresponds to Experiment 2 in Section [5] 
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(b) c = 100, 000. 

Figure 5: Illustration of Theorem 1321 Corresponds to Experiment 3 in Section [5] 



42 



